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Abstract 

We consider an arbitrary supersymmetric U(l) gauge factor at the TeV scale, 
under which the two Higgs superfields H± 2 of the standard model are nontrivial. We 
assume that there is a singlet superfield S such that H1H2S is an allowed term in 
the superpotential. We discuss first the generic consequences of this hypothesis on the 
structure of the two-doublet Higgs sector at the electroweak energy scale, as well as 
Z — Z' mixing and the neutralino sector. We then assume the existence of a grand 
unified symmetry and universal soft supersymmetry breaking terms at that scale. We 
further assume that the additional U(l) is broken radiatively via a superpotential term 
of the form hh c S, where h and h c are exotic color-triplet fields which appear in Eq 
models. We show that the U(l) breaking scale and the parameter tan/3 = V2/V1 are 
then both predicted as functions of the H1H2S coupling. 



1 Introduction 



If supersymmetry is broken at the TeV energy scale and the standard SU(3)c X SU(2)l X 
U(l)y gauge symmetry is not extended at that same scale, then the existence of supersymme- 
try above that scale protects the theory from nondecoupling contributions of physics above 
that scale. In other words, we get the Minimal Supersymmetric Standard Model (MSSM). 
However, if the gauge symmetry is extended also at the TeV energy scale and it breaks down 
to that of the standard model together with the supersymmetry, there will be in general new 
important phenomenological consequences at the TeV scale as well as the 100 GeV scale. 
With further simplifying assumptions, the parameters of the two scales may also be related. 

A particularly interesting extension of the MSSM is the inclusion of an extra U(l) factor 
at the TeV energy scale. The motivation for this could be theoretical. If the standard model 
is embedded in a larger symmetry group of rank greater than 4, such as 5*0(10) [rank 5] or 
Eq [rank 6], then an extra U(l) gauge factor is very possible at the TeV energy scale. This is 
especially true for the supersymmetric E 6 model |], || based on the E 8 x Eg heterotic string. 
In particular, if only flux loops are invoked to break E$ down to SU(3)c x SU(2)l x U(1)y, 
then a specific extra U(l) [conventionally known as U(l) v ] is obtained. Remarkably, U(l) v 
is also phenomenologically implicated || by the experimental = T(Z — > bb)/T(Z — > 
hadrons) excess. Another possible clue is the totality of neutrino-oscillation experiments 
(solar, atmospheric, and laboratory) which suggest that there are at least 4 neutrinos. This 
has been shown |4j] to have a natural explanation in terms of the E§ superstring model with 
a specific U(l) called Z7 (1)jv- 

In Sec. 2, we consider a generic extra supersymmetric U(l) gauge factor at the TeV 
energy scale with two doublet superfields H\ 2 and a singlet superfield S such that H1H2S 
is an allowed term in the superpotential. (Note that if S has a nonzero charge under the 
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aditional U(l), as is the case if the scalar component of S is to acquire a nonzero vacuum 
expectation value (VEV) so as to break this U(l), then above this breaking scale, no 11H1H2 
superpotential term exists.) We then derive its nondecoupling effects on the two-doublet 
Higgs sector at the 100 GeV scale. In Sec. 3, we specialize to a class of U(l) a models 
derivable from E$, of which the U(l) v and U(1)n models are special cases. In Sec. 4, we 
discuss how the new Z' mixes with the standard Z in the general case, and formulate the 
effects in terms of the oblique parameters €1,2,3 or S,T,U in the U(l) a models. We also 
discuss the generic neutralino sector. In Sec. 5, we show how supersymmetric scalar masses 
are affected by the extra D-terms from U(l) a . Combining this with the results of Sec. 2 and 
Sec. 3, and assuming universal soft supersymmetry breaking terms at the grand- unification 
scale, we show that there is a relationship between the U(l) a vacuum expectation value and 
the well-known parameter tan/3 = Vijv\ used in the MSSM. Finally in Sec. 6, we have some 
concluding remarks. 

2 Tree-Level Nondecoupling at the 100 GeV Scale 

As the U(l) gauge factor is broken together with the supersymmetry at the TeV scale, 
the resulting heavy scalar particles have nondecoupling contributions to the interactions 
of the light scalar particles. Consequently, the two-doublet Higgs structure is of a more 
general form than that of the Minimal Supersymmetric Standard Model (MSSM). Previous 
specific examples have been given. || |7], [8| Here we present the most general analysis. We 
denote the scalar components of Hi, H2, and S as $1, $2, and x respectively. Under 
SU(3) C x SU(2) L x U(1) Y x U(1) x , we then have 




(1) 
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$ 2 = | ^ J ~(l,2,i;-l + a), (2) 
X = x°~(l,l,0;l), (3) 

where each last entry is the arbitrary assignment of that scalar multiplet under the extra 
U(l)x with coupling g x , assuming of course that the superpotential has the term fH 1 H 2 S. 
The corresponding scalar potential contains thus 

V F = / 2 [($ t 1 $ 2 )($t$ 1 ) + ($t $1 + <!>l$ 2 )x X }, (4) 
and from the gauge interactions, 

Vd = ^|[(^>l<^i) 2 + (<^^ 2 ) 2 + 2(<^i<^ 1 )(^+<^ 2 ) - 4(<^i^ 2 )(^+<^ 1 )] 

o 

+ I^[-$t$ 1 + $t$ 2 ]2 

+ ^gli-a^-il-a^l^ + xx] 2 - (5) 

Let (x) = u, then \^2Rex is a physical scalar boson with m 2 = 2g 2 u 2 , and the (<&\<&i)\/2Rex 
coupling is v / 2-u(/ 2 — g 2 a). Hence the effective ($|$i) 2 coupling Ai is given by 

Ai = \(9i+9i)+ 9 y- 2{f2 ;f a)2 



X 
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-(gl + gl) + 2af-L (6) 



Similarly 



A 2 = i(^ + ^) + 2(l-a)/ 2 -^, (7) 



A3 = -\ 9 f + \ 9 i + f-^ (8) 

A 4 = —gl + f, (9) 

where the two-doublet Higgs potential has the generic form 

V = m^l^x + m^^s + m^^^l^s + ^^O + iA^^l^!) 2 

+^A 2 ($^ 2 ) 2 + V^iX^) + A 4 ($ t 1 $ 2 )($^ 1 ). (10) 



From Eqs. (6) to (9), it is clear that the MSSM is recovered in the limit / = 0. [Note that 
m i2 only after U(l) symmetry breaking and it would be proportional to / if universal 
soft supersymmetry breaking is assumed.] Let (0$ 2 ) = v h 2 -> tan/5 = v 2 /vi, and v 2 = vf + v%, 
then this V has an upper bound on the lighter of the two neutral scalar bosons given by 

(m 2 h ) max = 2v 2 [X 1 cos 4 (3 + A 2 sin 4 /3 + 2(A 3 + A 4 ) sin 2 /3 cos 2 /3] + e, (11) 

where we have added the radiative correction || due to the t quark and its supersymmetric 
scalar partners, i.e. 



3 9lmf ( fif 
S^M 2 , [ m 2 

We note also that the soft supersymmetry-breaking term fAf&\& 2 X + h.c. (from which 
we obtain m\ 2 = fAju) also contributes to A4 and generates some additional quartic scalar 
couplings. However, we assume here that fAf/g 2 u is small, because we are mainly interested 
in the case where the elctroweak Higgs sector has two relatively light doublets and not just 
one light doublet. Using Eqs. (6) to (9), we obtain 

P 



(m z h ) max = M l z cos z 2f3 + t + 
where 



A- J - 

9x 



(13) 



A= ^ + (2a-l)cos2/3-icos 2 2/?. (14) 

If A > 0, then the MSSM bound can be exceeded. However, f 2 is still constrained from the 
requirement that V be bounded from below. We note here that although Vp of Eq. (4) and 
Vd of Eq. (5) are nonnegative for any value of /, V of Eq. (10) is not automatically bounded 
from below. This simply means that if / is too large, the minimum of the original potential 
only breaks the extra U(l) but not the electroweak gauge symmetry. Given g x and a, we can 
vary cos 2(3 and / to find the largest numerical value of m^. We show in Fig. 1 this upper 
bound on as a function of g 2 for several specific values of a. The value ao is chosen in 
the top curve to maximize for a given value of g 2 . This upper bound increases as g 2 



increases. However, it is reasonable to assume that g x cannot be too large. In fact, in the 
specific models to be discussed in the next section, g 2 x < 0.16. As shown in Fig. 1, even for 
g 2 x = 0.5, the upper bound is only about 190 GeV. 

It should be mentioned that an upper bound on rrih has been previously obtained fTTfl 
assuming that there is no extra U(l) at the supersymmetry-breaking scale. However, the 
same proof also goes through with an extra U(l). We improve on Ref. [10] in this case by 
computing exactly how the off-diagonal nondecoupling terms affect the upper bound on rrih, 
resulting in Fig. 1 as shown. If fAf/g x u is not small as we have assumed, then the reduction 
to V of Eq. (10) is not valid. In this case, the electroweak Higgs sector consists of effectively 
only one Higgs doublet with rrih bounded by a function which involves also At .fill 



3 U(l) Gauge Factor from E(6) 

As already mentioned in the Introduction, an extra supersymmetric U(l) gauge factor at 
the TeV scale is a very viable possibility from the spontaneous breakdown of E§. Consider 
the following sequential reduction: 

E 6 -> 5O(10) [xf/(l) v> ] (15) 
5O(10) SU(5) [xC/(l) x ] (16) 
3U(5) -> SU(3) C x SU(2) L [xU(l) Y ]. (17) 

At each step, a U(l) gauge factor may or may not appear, depending on the details of the 
symmetry breaking. Assuming that a single extra U(l) survives down to the TeV energy 
scale, then it is generally given by a linear combination of U(l)^ and U(l) x which we will 
call U{l) a . 

Under the maximal subgroup SU(3)c X SU (3) l x SU(3)r, the fundamental representation 



6 



of Eq is given by 

27 = (3, 3, 1) + (3*, 1, 3*) + (1, 3*, 3). (18) 
Under the subgroup SU(5) x U(l)^ x U(l) x , we then have 

27 = (10; 1,-1) l(u,d),u c ,e c ] 

+ (5*; 1,3) [d c ,(v e ,e)\ 

+ (1; 1, -5) [N] 

+ (5; -2, 2) [h,(E c ,N c E )\ 

+ (5*; -2, -2) [h c ,(v E ,E)\ 

+ (i;4,o) [S], (19) 

where the U(l) charges refer to 2y / 6<3,/> an d 2-\/T0(5 x . Note that the known quarks and 
leptons are contained in (10; 1, —1) and (5*; 1,3), and the two Higgs scalar doublets are 
represented by [ve-.E') and (E C ,N E ). Let 

Qa = Q^p cos a — Q x sina, (20) 



then the rj- mo del |T], is obtained with tana = y3/5 and we have 

27 = (10;2) + (5*;-l) + (l;5) 

+ (5;-4) + (5*;-l) + (l;5), (21) 

where 2\/l5Q r] is denoted; and the iV-model[|[] is obtained with tana = —l/y/TE resulting 



in 



27 = (10;l) + (5*;2) + (l;0) 

+ (5; -2) + (5*; -3) + (1; 5), (22) 

where 2\/T0Qn is denoted. This model is so called because the superfield TV has Qn = 0. It 
allows S to be a naturally light singlet neutrino and is ideally suited to explain the totality 
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of all neutrino-oscillation experiments, i.e. solar|]12[, atmospheric [|l3f , and laboratory jnj. It 
is also a natural consequence of an alternative SO (10) decomposition ||T5|j of E$, i.e. 

16 = [(u,d),u c ,e c ;h c ,(v E ,E);S), (23) 
10 = [h,(E c ,N c E );d c ,(is e ,e)\, (24) 
1 = [N], (25) 
which differs from the conventional assignment by how the SU(5) multiplets are embedded. 

Identifying $i, $2, and x with the scalar components of (ue, E), (E c , N e ), and S of which 
we can choose one copy of each via a discrete symmetry 0] to be the ones with VEVs, we see 
that the general anaylsis of the previous section is applicable for this class of U(l)-extended 
models. (Of course, more than one copy of (u E ,E), (E c ,N e ), or S could have VEVs but 
that would lead to a much less constrained scenario.) Assuming that U(l) a is normalized in 
the same way as U(l)y, we find it to be a very good approximation || to have g 2 a = (5/3)gf. 
We then obtain for the ?7-model, 

9 2 x = 7^9 2 n a = \ ( 26 ) 



and for the iV-model, 



whereas in the exotic left-right model|| |i"4fl , 

2 _ (gl + gl)(l-siT?e w ) 2 _ 2 
9 '~ 4(1 -2 sin 2 6 W ) ' a " tan ^- ( 28 ) 

These three specific points have been singled out in Fig. 1. Furthermore, when we take the 
squark masses to be about 1 TeV we find the largest numerical value of in the U(l) a 
models to be about 142 GeV , as compared to 128 GeV in the MSSM, and it is achieved 
with 



2 A /3/5cos2/3 

ta " a = - 3-00^2/3 ■ (29) 



which is possible in the 77-model, i.e. tana = v/3/5 and cos 2(3 = — 1. 



4 Z - Z' and Neutralino Sectors 



The part of the Lagrangian containing the interaction of $1^ and \ with the vector gauge 
bosons Ai(i = 1,2,3), B, and Z' belonging to the gauge factors SU(2) L , U(1) Y , and U(l) x 
respectively is given by 

+ - f^A? - ^B" + ig x (l - a)Z'^2\ 2 



+ W-ig x Z>») X \\ 



(30) 



where are the usual 2x2 Pauli matrices. With the definition Z = ((72^3 — 9iB)/gz, where 



= ygt + gi, the mass-squared matrix spanning Z and Z' is given by 



Ml z , 



9z9x[-avf + (1 - a)u|] 



»zffx[-aui + (1 " 2^[ u 2 + + (1 - a) 2 *; 2 ] 



Let the mass eigenstates of the Z — Z' system be 



(31) 



Z\ — Z cos 9 + Z' sin Z 2 = — Z sin + Z' cos 6 1 , 



(32) 



then the experimentally observed neutral gauge boson is identified in this model as Zi, with 
mass given by 



and 



M 2 Zl =M* z c \g\J 



9z 



1 - (sin 2 (3 - a) 



u 



9 ~ -^(sin 2 /3-a)- 



(33) 



(34) 



2g x w 

Note that Z 2 has essentially the same mass as the physical scalar boson \[2Rex discussed 
earlier. 

So far, our discussion of the Z — Z' sector is completely general. However, in order to 
make contact with experiment, we have to specify how Z' interacts with the known quarks 



and leptons. In the class of U(l) a models from Eq, all such couplings are determined. In 
particular, we have 



2 1 I i 

9x = y -g a cosa, a = - I 1 



tan a 



(35) 



Using the leptonic widths and forward-backward asymmetries of Z\ decay, the deviations 
from the standard model are conveniently parametrized ||1 5 ] : 



£l 



£2 = -(3 — v 15 tan a) 



3 

- tana 
5 



- ~ aT, 
u z 



1 



sin' /3 — - 1 



f.3 



/3 1 

1 — 3W- tana H 5 

V 5 2 sin 2 9 W 



- tana 
5 



3 

- tana 
5 



aU 



u- 



4 sin 2 9 W ' 



(36) 
(37) 



x 



sin 2 /5- I I 1 



3 

- tana 
5 



aS 



(38) 



4 sin 2 W ' 

Since the experimental errors on these quantities are fractions of a percent, u ~ TeV is 
allowed. 

In the MSSM, there are four neutralinos (two gauge fermions and two Higgs fermions) 
which mix in a well-known 4x4 mass matrixfTT||. Here we have six neutralinos: the gauginos 
of U(l)y and the third component of SU(2)l, the Higgsinos of 4>1 and 0°> the U(l)x gaugino 
and the x Higgsino. The corresponding mass matrix is then given by 



Mi 





-<7i«i/V2 


giv 2 /V2 











M 2 


g 2 vx/\/2 


-g 2 v 2 / \p2 










g 2 vi/V2 





fu 


-g x aviy/2 


fV2 




-g 2 v 2 / V2 


fu 





-g x {l - a)v 2 V2 


M 








-g x aviV2 


-g x (l - a)v 2 \/2 


M x 


g x uV2 








fV2 


M 


g x u\f2 






(39) 

where M l 2 . 2 are allowed U(l) and SU(2) gauge-invariant Majorana mass terms which break 
the supersymmetry softly. Note that without the last two rows and columns, the above 
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matrix does reduce to that of the MSSM if fu is identified with — //. However, the \i 
parameter in the MSSM is unconstrained, whereas here fu is bounded and / itself appears 
in the Higgs potential. 

Since g x u should be of order TeV, the neutralino mass matrix M.^ reduces to either a 

4 x 4 or 2 x 2 matrix, depending on whether / is much less than g x or not. In the former 
case, it reduces to that of the MSSM but with the stipulation that the \i parameter must be 
small, i.e. of order 100 GeV. This means that the two gauginos mix significantly with the 
two Higgsinos and the lightest supersymmetric particle (LSP) is likely to have nonnegligible 
components from all four states. In the latter case, the effective 2x2 mass matrix becomes 

Mi + g\v x v 2 j 'fu -gig&xVi/ fu ^ 
-gi92Viv 2 /fu M 2 + g\v x v 2 l 'fu J 

Since V\V 2 ju is small, the mass eigenstates of M.'^ are approximately the gauginos B and 
W3, with masses M\ and M 2 respectively. In supergravity models with uniform gaugino 
masses at the GUT breaking scale, 

Mx = ^\M 2 ~ 0.5M 2 , (41) 

hence B would be the LSP, which makes it a good candidate for cold dark matter. 

5 Supersymmetric Scalar Masses 

The spontaneous breaking of the additional U(l) gauge factor at the TeV scale is not possible 
without also breaking the supersymmetry. 0] As a reasonable and predictive procedure, we 
will adopt the common hypothesis that soft supersymmetry breaking operators appear at 
the grand-unification (GUT) scale as the result of a hidden sector which is linked to the 
observable sector only through gravity. Moreover these terms will be assumed to be universal, 
i.e. of the same magnitude for all fields. 
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M' M = 



Consider now the masses of the supersymmetric scalar partners of the quarks and leptons: 

m\ = ml + m\ + m 2 F + m 2 D , (42) 

where mo is a universal soft supersymmetry breaking mass at the GUT scale, m\ is a 
correction generated by the renormalization-group equations running from the GUT scale 
down to the TeV scale, mp is the explicit mass of the fermion partner, and m 2 D is a term 
induced by gauge symmetry breaking with rank reduction and can be expressed in terms 
of the gauge-boson masses. In the MSSM, m 2 D is of order Mf and does not change m B 
significantly. In the f/(l) Q -extended model, m 2 D is of order Mf, and will affect m B in a 
nontrivial way. For example, in the case of ordinary quarks and leptons, 



Am 2 D (10;l,-l) = 




Am^(5*;l,3) = ^Mf, |l - 3^ jj tana 

This would have important consequences on the experimental search of supersymmetric 
particles. In fact, if mp is not too large, it is possible for the exotic scalars (which may 
be interpreted as leptoquarks depending on their Yukawa couplings) to be lighter than the 
usual scalar quarks and leptons. We have already discussed this issue in Ref. [18]. 

Assuming Eq. (42), we first consider the spontaneous breaking of U(l) a , i.e. (x) = u, 
which requires m 2 x to be negative. This may be achieved by considering the superpotential 

W = fH^S + f'hh c S + \ t H 2 Q 3 t c , (45) 

(where we have omitted the rest of the MSSM Yukawa couplings) together with the trilinear 
soft supersymmetry-breaking terms 

V soft = fA f Q\Q 2X + f'A f ,hh c X + \tA t $ 2 Q 3 t c (46) 

along with the soft supersymmetry-breaking scalar masses. Starting with a wide range of 
given values of m , the universal gaugino mass m^j, and the universal trilinear massive 
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parameter A at the GUT scale, we find that m 2 x does indeed turn negative near the TeV 
energy scale for many typical values of / and /'. An example of this is given in Fig. 2. 
The evolution of m 2 is mostly driven by /', but / also contributes primarily through its 
direct effect on Ay. From the negative value of m 2 at the TeV scale, we then obtain the 
predicted mass of Z', i.e. Mz> = (—2m 2 ,) 1 / 2 , which is also the mass of the physical scalar 
boson \/2Rex- However, as we will discuss shortly, the mass of the Z' so obtained must also 
be consistent with the desired electroweak symmetry breaking conditions. 

We assume that the top quark's pole mass is 175 GeV, and that at 1 TeV a s = 0.1 
which corresponds to a s (M z ) 0.12. We will also assume that at the TeV scale and above, 
the particle content of the model is that of three complete 27's of E e and some additional 
field content so as to achieve gauge coupling unification. The additional field content could 
be near the unification scale and hence provide threshold corrections that allow the gauge 
couplings to unify, perhaps even at the string compactification scale. Another possibility 
is to add an anomaly-free pair of SU (2) L doublet fields so as to mimic gauge coupling 
unification in the MSSM. Such an example is discussed for the a = N model of Ref. [8]. 
This model has the same unification scale as is possible in the MSSM. In calculating the 
gauge-coupling beta functions, we will in fact assume the field content of that model, but 
the choice of additional matter fields or threshold corrections to bring about gauge coupling 
unification has no significant effect on our calculation. The fact that such models have three 
complete 27's has the noteworthy implication that the gauge coupling at the unification 
scale is approximately the strong coupling. The reason is that with three copies of h and h c , 
the beta function for a s is zero in one loop above the TeV scale. Similarly, the gluino mass 
also does not evolve in this approximation. 

Defining 

V = 8tt 2 — — (47) 
a In// 
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(where /x is the scale), the relevant renormalization group equations are: 

Pin A 2 = -E^ 2 + 6A t 2 + / 2 , 

i 

Vlnf = -£^ 2 + 3A 2 + 4/ 2 + 3/' 2 , 

i 

Vlnf 2 = -£cF\ 2 + 3/ 2 + 5f 2 , 

i 

for the Yukawa couplings, 

VA t = £ c f ] 9iMi + 6A t 2 A t + / 2 A f , 

i 

VA f = J2^ ) 9 2 M l + 3X 2 t A t + 4fA f + 3f 2 A f , 

i 

VA f , = Y^g^M. + SfAf + BfAf,, 

i 

for the trilinear scalar parameters Ai, and 



Vm 2 s = 


-E^V + 2/ 2 X / + 3f 2 X 

i 


Vml = 


i 


Vml c = 


~Z^ c i 9i+ J x f , 

i 


Vm 2 ^ = 


-E^ 2 + / 2 ^, 


Vm% 2 = 


-E^ 2) ft 2 + 3A^ + / 2 ^/ 

i 


= 


-E4% 2 + a 2 ^, 

i 


Vrn^c = 


-E^ 2 + a 2 x,, 



where we have defined 

Xt = m 2 Qs + m 2 t c + m| 2 + A 2 t , 

^/ = m 5 + m li + m % 2 + ^/ ' 

X r = m 2 s + m 2 h + m 2 h c + A), , 
14 



and the coefficients cj field ' ) have the obvious values. Further, the gaugino mass M { scales the 
same as aij. These equations are modified in an obvious manner if tan/3 is large enough 
that At and A T cannot be ignored or if there are more than one sizeable coupling serving 
the purpose of /', which is certainly possible since we have three copies of h and h c in these 
models. 

A very important outcome of Eq. (42) is that the U(l) a and electroweak symmetry 
breakings are related. To see this, go back to the two-doublet Higgs potential V of Eq. (10). 
Using Eqs. (6) to (9) and Eq. (35), we can express the parameters m 2 2 , ml, and m\ in terms 
of the mass of the pseudoscalar boson m^, and tan/3. 



m \i = ~ m A sm P cos Pi 
ml = m 2 A sin 2 (3 - ^Mf cos 2/3 



2f 



2 sin 2 /3 + I 1 - J - tan a | cos 2 /3 % 

V 5 / 2 cos 2 a g z a 



(64) 



(65) 



m\ = m\ cos 2 (3 + -M| cos 2/3 



2/ 2 o 
9z 



2 cos 2 (3 + [ 1 + J ^ tan a | sin 2 (3 - - — %- tt 

V 5 / 2 cos 2 a g l a 



(66) 



On the other hand, using Eq. (42), we have 



™12 = f A f U i 



m 1 = m + m R1 + f u - - \1 - ^-tmin ] M z z ,. 
ml = mo + m| 2 + /V-| fl + W-tano;] M|,, 



(67) 
(68) 

(69) 



where m 2 Rl and m 2 R2 differ in that \ t (the Yukawa coupling of $2 to the t quark) contributes to 
the latter but not to the former. Both depend on m , m^, A , the various gauge couplings 
gi, as well as / and /'. Matching Eqs. (64) to (66) with Eqs. (67) to (69) allows us to 
determine u and tan/3 for a given set of parameters at the grand-unification scale. 
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We will now briefly discuss our method for finding u and tan (3 for a given set of univer- 
sal soft supersymmetry-breaking parameters mj,mo, Aq at the GUT scale and the Yukawa 
coupling /, when such a solution exists. First, we guess a value for tan/? so as to choose 
a value for X t . We then form a table [M z >,m 2 R1 ,m 2 R2 , Af\(f ) for many very closely spaced 
values of /' extending up to where f'(Ma) reaches its perturbative limit. By "closely spaced 
values of /'," we mean that between two consecutive entries in the table, none of the four 
parameters differs by more than one percent. Second, we guess a value for M z > which lies 
within the range in the table, so as to choose m 2 R1 , m 2 R2 , Af from the entry of the table which 
has M Z i closest to this value. Third, we equate the right-hand sides of Eq. (65) + Eq. (66) 
and Eq. (68) + Eq. (69) to solve for m 2 A as a linear function of u 2 and cos 2 /5. Fourth, 
using the previous result for m 2 A we equate the right-hand sides of Eq. (65) — Eq. (66) and 
Eq. (68) — Eq. (69) to solve for u 2 as a function of cos 2 (3 of the form of a linear function 
divided by another linear function. Fifth, using the expressions from the previous two steps 
we equate the right-hand sides of Eq. (64) with that of Eq. (67) and solve numerically for 
cos 2 f3, and hence tan f3, by first searching for a root close to the value corresponding to our 
original guess for tan/5. In doing this fifth step, one needs to choose fu > or fu < 
analagous to /i > or /j, < in the MSSM, and then check that the solution is consistent 
with m\ = —fAfu/ sin (3 cos (3 > 0. In fact, taking all Yukawa couplings and tan/5 to be 
positive as well as our convention for the trilinear coupling parameters Ai, solutions exist 
only for u < 0. Next, if a solution to these steps has been found, we start the entire cycle 
over using the values for tan/5 and M Z i just calculated as the new "guessed" values. This 
iteration is continued until the predicted tan (3 and M z > become fixed to a reasonable accu- 
racy (we demand about five-percent accuracy). This process can be speeded up by adding 
a sixth step to the cycle which repeats the third through fifth steps until the prediction for 
tan (3 and M z > become fixed for the table found in the second step of the cycle. 
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Before we discuss our results, we remind the reader that for the case that Af could be 
small / has a maximum possible value that comes from requiring that the Higgs potential 
be bounded from below and which depends on the additional U(l). We plot this maximum 
value / max as a function of a (see Eq. (20)) in Fig. 3. In particular, the r\— model requires / 
to be less than about 0.35, whereas for a — 0, / could be as great as 0.46. Note that as |a| 
approaces n/2, / max approaches 0. From Fig. 2(a), one can see that if / is small enough so 
that a = i] is allowed, then f(Ma) will always be perturbatively small for a perturbatively 
valued /'. In our examples, we will only be interested in values of / < 0.35. 

In Fig. 4, we show the predicted values of tan (3 and M z < as a function of a for / = 0.345 
and rrig = 200 GeV, A = 650 GeV and m = 650 GeV. In accordance with Fig. 3, we 
are only interested in showing |a| less than about 0.7. We have also plotted the magnitude 
\u\ of the VEV of the singlet Higgs boson x and the mass f'\u\ of the exotic fermion h(h c ). 
In Fig. 5, we show the similar situation for / = 0.345 and m g = 300 GeV, A = 950 
GeV and m = 950 GeV. These two figures are quite similar except that the mass scale in 
Fig. 5(b) has been pushed up relative to that shown in Fig. 4(b). These choices of soft 
supersymmetry-breaking parameters are fairly typical in that generally we need m Q to be at 
least twice as great as m g to find a solution. Further, if we want to have a solution for all a 
less than some value, A must be positive and of order m . 

In Figs. 6-9, we illustrate the effects of varying the parameters f,m g ,A Q ,m Q for a fixed 
value of a — rj. We look at the solutions for tan (3 and M z > (as well as \u\ and f'\u\) when 
the four input parameters are varied one at a time around the point / = 0.345 and m g = 250 
GeV, A = 650 GeV and m = 650 GeV. Note from Fig. 6 that with decreasing /, tan/5 and 
M Z i both increase. We do not extend / above 0.345 so as to avoid the upper bound coming 
from Fig. 3. We find also that we cannot decrease / much below 0.32 for this example and 
still have a solution for the electroweak breaking. To use smaller values of /, one would have 
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to increase the scale of the soft supersymmetry-breaking parameters. In Fig. 7, we look at 
the effect of varying Aq. The range of Aq examined is restricted because any extension in 
either direction would require values of M z > larger than can be reached via the f'hh c S term 
with /' within the perturbative regime. In Fig. 8, we vary m . Note that with increasing m , 
the predicted tan (3 increases significantly and M Z r decreases. In this example, increasing ijiq 
beyond 1200 GeV would predict an M z > less than 500 GeV and a tan/5 greater than 10. The 
lower limit of 500 GeV for m used here is due to the same reason as just given for the range 
of Aq plotted in the previous figure. In Fig. 9, we show the effect of varying the gluino mass 
which is also here the GUT scale universal gaugino mass. With increasing gluino mass, tan (3 
decreases while M z > increases. The upper limit of 350 GeV used here for the gluino mass 
again corresponds to about the size of that parameter for this example where increasing it 
anymore would require values of \u\ larger than can be reached perturbatively through the 
renormalization-group equations. We find the general trends of Fig. 6-9 to be typical of 
other choices of paramter values where consistent solutions exist. 

If m is demanded to be less than about 1 TeV, then in general tan (3 < 10, where the b and 
r Yukawa couplings are small enough not to contribute significantly to the renormalization- 
group equations. It is interesting to note that in contrast to the MSSM, where m\ — m 2 , = 
—m 2 R2 (\ t ) = — (m 2 A +m 2 z ) cos 2/3, solutions with tan f3 < 1 in principle are possible here due to 
the TeV scale D-terms. However, to have such a solution in practice with m| pole ' ) m 175 GeV 
means having \t(m t ) greater than its fixed-point value of about 1.22 with «g = a s (l TeV) ~ 
0.1 where the gauge couplings run according to the additional exotic field content as we have 
chosen. 

If fAf/g^u, where g 2 = (2/3)gj cos 2 a is not small, then Eqs. (64) to (66) have additional 
contributions, but they are always suppressed by v 2 /u 2 relative to m\ 2 = fAfU, hence our 
numerical results on tan/5 and M z >, etc. do not change appreciably. The corrections are 
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only important if the masses and splittings of the two Higgs doublets are considered. 

6 Conclusions 

We have shown in this paper that there are many interesting and important phenomenological 
consequences if we assume the existence of a supersymmetric U(l) gauge factor at the TeV 
energy scale. We assume that there is a Higgs superfield S which is a singlet under the 
standard gauge group but which transforms nontrivially under this extra U(l) so that it 
may break the latter spontaneously without breaking the former. We assume also that 
HiH 2 S is an allowed term in the superpotential. We then analyze the most general form of 
the Higgs potential and derive an upper limit on the lighter of the two neutral scalar Higgs 
bosons of the two-doublet Higgs sector as shown in Fig. 1. This generalizes the well-known 
case of the Minimal Supersymmetric Standard Model (MSSM). 

We then specialize to the case where this extra U(l) is derivable from a E 6 model with the 
particle content given by its fundamental 27 representation. We discuss the effect on Z — Z' 
mixing and the oblique parameters 61,2,3, as well as the extended neutralino mass matrix. We 
then work out in detail the consequences for supersymmetric scalar masses. We note that 
the mere existence of a spontaneously broken U(l) gauge factor at the TeV scale implies 
new important corrections to these masses through the so-called D-terms which are now 
dominated by Mf, instead of just Mf in the MSSM. This changes the entire supersymmetric 
scalar particle spectrum and should not be overlooked in future particle searches. 

Assuming universal soft supersymmetry-breaking terms at the grand-unification (GUT) 
scale, we match the electroweak breaking parameters with the corresponding ones from the 
U(l) breaking. Specifically the values of mf, m%, and m\ 2 in the well-known two-doublet 
Higgs potential are constrained as shown by Eqs. (64) to (69). We then obtain consistent 
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numerical solutions to these constraints and demonstrate how the U(l) breaking scale and 
the parameter tan/3 = Vi/vi are related through the HiH 2 S coupling. Our results are 
presented in Figs. 2 to 9. 

During the final stage of completing this manuscript, we became aware of Ref. |TD[ , which 
also discusses electroweak symmetry breaking with an additional supersymmetric U(l) gauge 
factor, but the emphasis there is on the case /' = 0. The case /' ^ is also discussed there, 
but the conclusion is that whereas the breaking of the additional U(l) radiatively via the 
term f'hh c S, already noted in Ref. ||, can be achieved with universal soft supersymmetry- 
breaking terms at the GUT scale, it does not work in the large trilinear coupling scenario. 
Our approach is essentially orthogonal. We concentrate on solutions where the U(l) scale is 
much larger than the electroweak scale. With the two scales being intimately related through 
the matching of Eqs. (64) to (66) with Eqs. (67) to (69), it is in fact highly nontrivial to find 
solutions which are consistent with this matching even with an arbitrary /'. We note also 
that our examples are models with complete Eq particle content and in our approximation, 
the Yukawa coupling / is bounded as shown in Fig. 3. In the more general case, the bound 
on / increases as the trilinear coupling increases. 
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Figure Captions 



Fig. 1. The upper bound on the lighter Higgs mass rrih as a function of g\ for various values 
of a. In all cases, we find the allowed value of / = f that maximizes m h . In the top curve, 
we find the pair f — fo and a = a that maximizes rrih whereas the value of a is held fixed 
as labeled for the other curves. The points corresponding to the 77, N, and exotic left-right 
models, described in Section 3, are marked by arrows. 

Fig. 2(a). The parameter / at 1 TeV as a function of f G = f(M G ) for models with an extra 
U(l) originating from E e . In descending order, the curves represent f' G = f'(M G ) = 0.5, 1.0, 
2.0, 3.0. 

Fig. 2(b). The mass M z > as a function of fa with the same values m g = 250 GeV, A = 650 
GeV and m = 650 GeV for different curves with the values of f' G as in 2(a). 

Fig. 3. The maximum value of / = / max for which the Higgs potential is bounded from below 
as a function of a, defined in Eq. (20). 

Fig. 4(a). tan/3 as a function of a for / = 0.345 and m g = 200 GeV, A = 650 GeV and 
m = 650 GeV. 

Fig. 4(b). Mz' (solid line), \u\ (short-dashed line) and f'\u\ (long-dashed line) as a function 
of a for the same values of input parameters as in 4(a). 

Fig. 5(a). tan/3 as a function of a for / = 0.345 and m~ g = 300 GeV, A = 950 GeV and 
mo = 950 GeV. 

Fig. 5(b). Mz' (solid line), |w| (short-dashed line) and f'\u\ (long-dashed line) as a function 
of a for the same values of input parameters as in 5(a). 

Fig. 6(a). tan/3 as a function of / for a = rj and m~ g = 250 GeV, A = 650 GeV and 
m = 650 GeV. 
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Fig. 6(b). M z > (solid line), \u\ (short-dashed line) and f'\u\ (long-dashed line) as a function 
of / for the same values of input parameters as in 6(a). 

Fig. 7(a). tan/? as a function of A for a = rj and / = 0.345, m g = 250 GeV and m = 650 
GeV. 

Fig. 7(b). M z > (solid line), \u\ (short-dashed line) and f'\u\ (long-dashed line) as a function 
of Aq for the same values of input parameters as in 7(a). 

Fig. 8(a). tan/? as a function of m for a = rj and / = 0.345, m~ g = 250 GeV and A = 650 
GeV. 

Fig. 8(b). Mz> (solid line), \u\ (short-dashed line) and f'\u\ (long-dashed line) as a function 
of mo for the same values of input parameters as in 8(a). 

Fig. 9(a). tan (3 as a function of rrig for a = r\ and / = 0.345, A = 650 GeV and m = 650 
GeV. 

Fig. 9(b): Mz< (solid line), |w| (short-dashed line) and f'\u\ (long-dashed line) as a function 
of mo for the same values of input parameters as in 9(a). 
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